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The S=l/2 antiferromagnetic Heisenberg chain exhibits 
a magnetic field driven quantum critical point. We study 
the low frequency spin dynamics in copper pyrazine 
dinitrate (CuPzN), a realization of this model system 
of quantum magnetism, by means of 13 C-NMR spec- 
troscopy. Measurements of the nuclear spin-lattice relax- 
ation rate Tj - in the vicinity of the saturation field are 
compared with quantum Monte Carlo calculations of the 



dynamic structure factor. Both show a strong divergence 
of low energy excitations at temperatures in the quantum 
regime. The analysis of the anisotropic T 1 ~ 1 -rates and 
frequency shifts allows one to disentangle the contribu- 
tions from transverse and longitudinal spin fluctuations 
for a selective study and to determine the transfer of de- 
localized spin moments from copper to the neighboring 
nitrogen atoms. 
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1 Introduction Low dimensional quantum magnets, 
e.g. spin chains or ladders, are model systems for the 
study of quantum criticality 1 1 2 3]|4][5]|2]|2][EE) • ^ or 
the isotropic S=l/2 antiferromagnetic Heisenberg chain 
(AFHC), nuclear magnetic resonance (NMR) has been 
used to probe the slow spin fluctuations in the full magnetic 
phase diagram of the AFHC, i.e., from the low to the high 
field limit and at temperatures from the quantum regime 
(k B T < J) up to the classical regime (k B T > J)fT0l. 
The Hamiltonian of the AFHC in an external field reads: 



H = ^( JSi ■ S l+ i - gnB~B ■ Si), 



(1) 



where Sj are spin operators and J is the exchange energy. 
At B c = 2 J / (ff/is) the system exhibits a quantum critical 
phase transition, where a single Ising triplet as lowest ele- 
mentary excitation crosses the ground state and the system 
switches from complete polarization into a Luttinger liquid 
of deconfined spinons. The condensation of low energy ex- 
citations is often referred to as critical slowing down and 
gives rise to a divergence of slow spin fluctuations. 



In this paper we present a NMR study of the low frequency 
spin response for various magnetic fields close to B c = 14.9 
T IfTTI in the metalorganic AFHC copper pyrazine dinitrate 
(CuPzN). We discuss experimental data for fields from 2.0 
to 14.8 T and temperatures from 1 .6 to 30 K. The dynamics 
are probed by the nuclear spin-lattice relaxation rate T-f" 1 
and compared with quantum Monte Carlo (QMC) calcula- 
tions. We also study the anisotropy of contributions from 
transverse and longitudinal spin fluctuations to T^ 1 and 
relate them to the anisotropic NMR frequency shift in the 
low field regime. The results are used to experimentally se- 
lect the transverse spin fluctuations in the present work and 
to evaluate the spin density transfer from the copper ions 
to the neighboring nitrogen atoms in the pyrazine ring. 

2 Experimental The compound CuPzN, i.e., Cu(C4- 
H4N2)(NC<3)2 has been characterized by inelastic neutron 
scattering, muon-spin relaxation, magnetothermal trans- 
port, specific heat and magnetization measurements IfTTI 
[T2l[T3l[T4l . All of these studies confirm CuPzN to be one 
of the best realizations of the AFHC with a low coupling 
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Figure 1 Segments of two neighboring chains A and B 
in the unit cell of CuPzN. The 13 C atoms reside in the 
pyrazine rings which are twisted against each other by ~ 
80° in the b-c plane. 

constant J/fcs = 10.7 K. This low value of J allows the 
study of the magnetic QCP in standard laboratory mag- 
nets. The electronic S =1/2 spins are mainly located in the 
Cu 2+ 3d- orbitals, whereas the 13 C nuclei in the pyrazine 
ring were used as / = 1/2 NMR probes. The onsite cop- 
per nuclei exhibit an extremely short spin-correlation time 
that is below the dead time of a standard NMR spectrome- 
ter. The measurements were performed in an 8 and a 17 T 
superconducting magnet with a modified Bruker CXP200 
spectrometer and a home-built spectrometer, respectively. 
In both cases a standard inversion-recovery spin-echo pulse 
sequence was used. 

3 Nuclear Magnetic Resonance The nuclear spin- 
lattice relaxation rate measures the electronic spin 
fluctuations at the nuclear Larmor frequency oj n . For the 
AFHC it can be written as 031 
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E E \ A lM+ A lM\ 



x / < Sp(q,t)Sp(-q,0) > 



l dt (2) 
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J2[F±(q)S ± (q,u> n ) + F z (q)Sz(q,uJ n )}. (3) 



-(<?) = Ea= 



[ A lp(l) 



Al p {q)}<mAF z {q) 



Here, F± 

A xz(o) + Ay Z (q) are the geometrical form factors and 
S±(q,uj) and S z (q,u>) are the dynamical structure factors 
of the electronic spin system. A a p with a, /3 = x,y,z are 
the components of the hyperfine coupling tensor A(q). 

3.1 Anisotropy of the NMR shift d In CuPzN both 
an isotropic hyperfine part A is0 (q) and an anisotropic 
dipolar part A dip (g) contribute to the coupling of the 13 C 
nuclei to the magnetic moments of the Cu 2+ electrons. 
The NMR shift S is proportional to the projection of the 
hyperfine coupling tensor A(o = 0) on the direction of 



the external field, i.e., the component A zz (0). Thus, for a 
rotation of the crystal in a fixed external field the shift can 
be extracted from: 



D 



(4) 



A(0)=D-((A dip (0) + A iso (0))-Xei- 

where Xei is the tensor of the electronic susceptibility, a 
is the diamagnetic tensor, and D is a rotation matrix. Fig. 
[2 a) shows the experimental shift 5 {(3) versus a simulation 
based on localized dipole moments for a rotation around 
the crystallographic a-axis, where f3 = 0° corresponds 
to an external field parallel to the o-axis, compare Fig. Q] 
The assumption that 100% of the spin magnetic moment 
resides on the Cu sites cannot describe the experimental 
data. A very good agreement is achieved for a spin moment 
transfer of 10% from copper to each of the neighboring ni- 
trogen atoms in the pyrazine ring and a misalignment of 5° 
between the actual rotation axis and the crystallographic 
a-axis. The transfer of delocalized spin moments was ob- 
served in similar systems lfl6l . For the resulting isotropic 
hyperfine part we find <5; so = 3600 ppm. The existence of 
two sets of spectral lines results from the relative twisting 
of pyrazine rings on two neighboring chains by ~ 80°. 
According to Fig.[T] those will be labeled as line set A and 
B in the following. 

3.2 Transverse dynamic structure factor Be- 
fore turning to the data, we present our method 
of calculation for the field and temperature-dependent 
transverse dynamic structure factor S±(q, u> n ). Switch- 
ing to imaginary time r the latter reads S±(q,r) = 
— J °° du>S±(q,ui)K(u),T), with a kernel K(ui,t) = 
e~ ru + e~^~ r ^ and (3 = 1/T. In real-space S±(q,T) 
can be calculated efficiently, using QMC. Following Ref. 

El 



Si,j (r) = ( E 



\n-m | 



" B n (n + l)(n- 



m)!m! 
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-(p)S-{p + m,)) 



(5) 



where S±(q,r) = £ a e i90 5 a ,o(r) /N and a, label lattice 
sites in a chain of length N. (■ ■ -)w refers to the Metropolis 
weight of an operator string of length n generated by the 
stochastic series expansion of the partition function Ifl8l 
[T9l , and p, m are positions in this string. 

Analytic continuation from imaginary times, i.e., 
S±(q,r), to real frequencies, i.e., S±(q,uj), is performed 
by the maximum entropy method (MaxEnt), minimizing 
the functional Q = \- 2 /2 - aa B20II211 . Here x refers 
to the covariance of the QMC data to the MaxEnt trial- 
spectrum S a ±(q, ui). Overrating is prevented by the en- 
tropy g = S a ±(q,u)hi{S a ±{q,uj)/m(uj)}. We have 
used a flat default model to(cj), matching the zeroth mo- 
ment of the trial spectrum. The optimal spectrum follows 
from the weighted average 



S±(q,u>)= / daP[a\S(q,T)]S a j_{q,uj) 



(6) 



Copyright line will be provided by the publisher 



pss header will be provided by the publisher 



3 




a) B = 2.0 T, T = 7 K. 
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b) B = 2.0 T, T = 7 K. 
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c) B = 2.0 T, T = 7 K. 



Figure 2 (a) NMR shift S versus rotation angle /3. As in (b) and (c), full (line set A) and dashed (line set B) lines 
correspond to calculations with a spin moment transfer of 10% in the simulation, gray lines correspond to a transfer of 0%. 
(b) Experimental T] -1 -rates versus f3 compared with a modulation of T-j - based on Eqs. (0 and (@}. (c) The modulation 
shown in (b) is disentangled into rates that stem from transverse (large values) and longitudinal (small values) fluctuations 
of the electronic moments. The arrow marks j3 = 50°, where the latter rate of set A goes to zero. 



with the probability distribution P[a\S(q,r)] adopted 
fromRef. (20). 

3.3 Anisotropy of spin fluctuations The ratio of 
longitudinal and transverse spin fluctuations contributing 
to the total spin-lattice relaxation rate at a given field can 
be evaluated from the anisotropy of T^ 1 l22l . Our previ- 
ous studies show that in this system the geometrical form 
factors for the case of 13 C-NMR are almost q-independent 
iflOl . Therefore, the spin-lattice relaxation rate only de- 
pends on S±(q), S z (q) and the orientation between the 
external field and the local hyperfine fields, which in turn 
are well known from the study of S (/?). Fig. 0(b) shows 
the experimental T^ 1 -rates versus a modulation of the 
relaxation rates based on equations (0, (01 and an almost 
isotropic dynamic structure factor calculated by QMC. 
The modulation was fitted to the experimental data by a 
single scaling factor determined by a least-squares fit. We 
find a good agreement for the comparison between our 
experimental data and the quantum Monte Carlo results. 
In this work we want to compare the experimentally and 
theoretically determined field- and temperature-dependent 

-rates stemming from transverse fluctuations only. 
Therefore, the contribution from longitudinal fluctuations 
S z (q,uj) was minimized by minimizing F z (q). Fig. 0(c) 
shows that for line set A the longitudinal fluctuation con- 
tribution to T^ 1 approaches zero for (3 = 50°. This line 
set and orientation was used for the subsequent measure- 
ments shown in Figs.0and0 Note that the transverse spin 
fluctuations also dominate the nuclear relaxation due to the 
large Fermi contact term in F±(q). 

3.4 1/Ti: field dependence In Fig. we compare 
the observed NMR rates with the QMC results versus 
magnetic field in the quantum regime fcgT <C J, with 
T = 1.6 K. The QMC data shown in Figs. and were 
calculated from the transverse dynamic structure factor 



S±(q,uj) l23l and scaled with a factor assigned at 2.0 T 
and high temperatures. The similarity between experiment 
and theory is remarkable. For both we find a pronounced 
maximum of T^ 1 (B) at B = 13.8 T shifting to lower 
fields with increasing temperature. To interpret these re- 
sults, we note that in the fully polarized state for B > B c , 
single magnons are exact eigenstates of Eq. (01 with a 
dispersion of 

E y (k) = J cos(fc) +gfi B B . (7) 

Ey (fc) displays a field-driven excitation gap of g^sB — J 
leading to an exponential decrease in T 1 _1 (_B) at fixed T 
and for B > B c . The rates calculated by QMC display a 
broader maximum than the measured data, but drop with 
the same slope for fields above 16 T. At B = B c the 




B[T] 

Figure 3 Field dependence of the nuclear spin-lattice re- 
laxation rate of 13 C close to the critical field. The expo- 
nential decrease of data reflects the linear opening of the 
spin gap with field. The inset highlights the maximum of 
T^ 1 {B) at B = B c and T = 1.6 K. All solid lines are a 
guide to the eye. 
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Figure 4 Temperature dependence of the nuclear spin- 
lattice relaxation rate of 13 C for different external fields 
in the vicinity of the quantum critical point. All QMC data 
lines are polynomial fits to analytic continuations of a 128 
site system. The theoretical errors range from 10% at low 
temperatures to 20% at high temperatures, estimated from 
different MaxEnt regularization parameters. 

dispersion touches the zero at k = ir/2 with a quadratic 
momentum dependence yielding a van-Hove type of crit- 
ical DOS. This leads to the maximum in T± , tending to 
diverge as T — > 0. For both, NMR experiment and QMC, 
the maximum in Fig. [3] occurs at B c 13.8 T, which is 
slightly less than the saturation field of B c — 14.9 T for 
the magnetization. This downshift is most likely a finite- 
temperature effect of excitations populating the gap. 

3.5 1/T X : temperature dependence In Fig. g] the 
temperature dependent rates for fields in the vicin- 
ity of the saturation field are shown. We find the strongest 
low temperature divergence at £? c =13.8 T. This is very sug- 
gestive of critical scattering as T — > 0. With increasing 
temperature the van-Hove singularity in the DOS at B c is 
smeared leading to the decrease in Tj~ . At B c ~ B = 
14.8 T the experimental data show no divergence but rather 
the opening of a spin gap. Possibly this field is slightly 
above the true saturation field for this orientation due to 
the sensitivity of the exact value of B c to the anisotropic g- 
factor. The QMC data for 12.8 and 13.8 T show a satisfying 
agreement with the experiment, with the same divergence 
at low temperatures. In contrast to the experimental data, 
at 14.8 T there is still a divergence as the quantum regime 
is approached. 

In the classical regime ksT J the rate is decreasing 
with increasing fields. This is indicative of an excitation 
spectrum dominated by spin-diffusion modes from q = 0. 

4 Conclusion We performed a complementary ex- 
perimental and theoretical study of the low-frequency spin 
dynamics in the isotropic S=l/2 antiferromagnetic Heisen- 



berg chain system copper pyrazine dinitrate. We find clear 
evidence for critical dynamics close to a field-induced 
QCP, as probed by the nuclear spin-lattice relaxation rate 
T\~ 1 . Experimental data and QMC calculations are in good 
agreement and show a pronounced maximum in Tf 1 in the 
vicinity of the saturation field. In this study, the anisotropic 
contribution from longitudinal spin fluctuations to T^ 1 
was minimized. From the analysis of the NMR frequency 
shift we conclude a spin moment transfer of 10% from 
copper to each of the neighboring nitrogen atoms in the 
pyrazine ring. 

Acknowledgements Part of this work was performed at the 
National High Magnetic Field Laboratory, supported by NSF Co- 
operative under Agreement No. DMR-0084173, by the State of 
Florida, by the DOE and the DFG under Grants No. KL1086/6- 
2 and No. KL 1086/8-1. One of us (W.B.) acknowledges partial 
support by the DFG through Grants No. BR 1084/4-1 and BR 
1084/6-1. He also thanks the KITP for hospitality, where this re- 
search was supported in part by the NSF under Grant No. PHY05- 
51164. 

References 

[1] B. C. Watson, V. N. Kotov, M. W. Meisel, D. W. Hall, G. E. 

Granroth, W. T. Montfrooij, S. E. Nagler, D. A Jensen, R. 

Backov, M. A. Petruska, G. E. Fanucci, and D. R. Talham, 

Phys. Rev. Lett. 86, 5168 (2001). 
[2] T. Lorenz, O. Heyer, M. Garst, F. Anfuso, A. Rosch, Ch. 

Riiegg, and K. Kramer, Phys. Rev. Lett. 100, 067208 (2008). 
[3] Z. Honda, H. Asakawa, and K. Katsumata, Phys. Rev. Lett. 

81, 2566 (1998). 
[4] V. S. Zapf, D. Zocco, B. R. Hansen, M. Jaime, N. Harrison, 

C. D. Batista, M. Kenzelmann, C. Niedermayer, A. Lacerda, 

and A. Paduan-Filho, Phys. Rev. Lett. 96, 077204 (2006) 
[5] S. A. Zvyagin, J. Wosnitza, C. D. Batista, M. Tsukamoto, N. 

Kawashima, J. Krzystek, V. S. Zapf, M. Jaime, N. F. Oliveira 

Jr., and A. Paduan-Filho, Phys. Rev. Lett. 98, 047205 (2007) 
[6] H. Manaka, I. Yamada, Z. Honda, H. Aruga Katori, and K. 

Katsumata, J. Phys. Soc. Jpn. 67, 3913 (1998). 
[7] V. O. Gaiiea, A. Zheludev, T. Masuda, H. Manaka, L.-P. 

Regnault, E. Ressouche, B. Grenier, J.-H. Chung, Y. Qiu, 

K. Habicht, K. Kiefer, and M. Boehm, Phys. Rev. Lett. 98, 

167202 (2007). 

[8] M. Takigawa, N. Motoyama, H. Eisaki, and S. Uchida, Phys. 

Rev. Lett. 76, 4612 (1996). 
[9] G. Chaboussant, P. A. Crowell, L. P. Levy, O. Piovesana, A. 

Madouri, and D. Mailly, Phys. Rev. B 55, 3046 (1997). 
[10] H. Kuhne, H.-H. Klauss, S. Grossjohann, W. Brenig, F. J. 

Litterst, A. P. Reyes, P. L. Kuhns, M. M. Turnbull, and C. P. 

Landee, Phys. Rev. B 80, 0451 10 (2009) 
[11] P. R. Hammar, M. B. Stone, D. H. Reich, C. Broholm, P. 

J. Gibson, M. M. Turnbull, C. P. Landee and M. Oshikawa, 

Phys. Rev. B 59, 1008 (1999). 
[12] T. Lancaster, S. J. Blundell, M. L. Brooks, P. J. Baker, F. L. 

Pratt, J. L. Manson, C. P. Landee and C. Baines, Phys. Rev. 

B 73, 020410(R) (2006). 
[13] M. B. Stone, D. H. Reich, C. Broholm, K. Lefmann, C. 

Rischel, C. P. Landee and M. M. Turnbull, Phys. Rev. Lett. 

91, 037205 (2003)] 



Copyright line will be provided by the publisher 



pss header will be provided by the publisher 



5 



[14] A. V. Sologubenko, K. Berggold, T. Lorenz, A. Rosch, E. 

Shimshoni, M. D. Phillips and M. M. Turnbull, Phys. Rev. 

Lett. 98, 107201 (2007). 
[15] T. Moriya, Prog. Theor. Phys. 16, 23 (1956). 
[16] A. U. B. Wolter, P. Wzietek, S. Sullow, F. J. Litterst, A. Ho- 

necker, W. Brenig, R. Feyerherm and H.-H. Klauss, Phys. 

Rev. Lett. 94, 057204 (2005). 
[17] A. W. Sandvik, J. Phys. A 25, 3667 (1992). 
[18] A. W. Sandvik, Phys. Rev. B 59, R14157 (1999). 
[19] O. F. Syljuasen and A. W. Sandvik, Phys. Rev. E 66, 046701 

(2002). 

[20] J. Skilling and R. K. Bryan, Mon. Not. R. Astron. Soc. 211, 
111 (1984). 

[21] M. Jarrell and J. E. Gubernatis, Phys. Rep. 269, 133 (1996). 
[22] L. J. Azevedo, A. Narath, P. M. Richards and Z. G. Soos, 

Phys. Rev. B 21, 2871 (1980). 
[23] S. Grossjohann and W. Brenig, Phys. Rev. B 79, 094409 

(2009) 



Copyright line will be provided by the publisher 



